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1 Introduction 



It has been suggested that occurrence of nonequilibrium chiral phase transi- 
tion in hadronic collisions would lead to the formation of disoriented chiral 
condensate (DCC), domains of coherent pion fields |1J. It then may offer a 
viable explanation of the Centauro events found in cosmic ray experiments 
which have anomalously large fluctuations of neutral to charged pion ratio . 

By a numerical simulation using the linear sigma model, Rajagopal and 
Wilczek || demonstrated that low-momentum components of pion fields are 
greatly amplified, indicating the formation of large coherent domains of pions. 
A key ingredient in their simulation is the "quench" initial condition which 
models a hypothesized rapid cooling of hot debris formed during the hadronic 
reactions. The formation of large domains of pion fields was confirmed by more 
realistic but still the sigma model simulation by Asakawa, Huang, and Wang 
0] in which the longitudinal expansion of hadronic blob is taken into account. 

There remain, however, questions regarding the interpretation of the sim- 
ulations. Among other things, the amplification of the low-momentum modes 
in the Rajagopal- Wilczek simulation lasts so long, much longer than l/m CT , 
the time scale of rolling down from the top of the Mexican hat potential. In 
previous works |S], £3J, we have conjectured that the parametric resonance mech- 
anism might be relevant for the DCC formation. In simple physical terms, the 
parametric resonance mechanism stands for a mechanism of amplification of 
pion field fluctuations due to cooperative effects coupled with oscillations of 
background sigma model fields. 

It has been known since long time ago that the parametric resonance mech- 
anism can be relevant for cosmology by providing reheating mechanism during 
inflation |7]]. It has been received renewed interests recently in relationship with 
various issues, such as the gravitino problem, in inflationary cosmology f8fl-|10|l. 
The parametric resonance mechanism for pion production was first discussed 
by Mrowczyhski and Miiller in close analogy to the out-of-equilibrium phase 
transition in cosmology [II]. We have stressed in our previous paper j5| that 
it may offer an efficient mechanism of pion production with potential possibil- 
ity of explaining long-lasting amplification of low momentum modes. We also 
noted that it would give an alternative picture of large isospin fluctuations, the 
picture quite different from the conventional one based on rolling down into the 
"wrong" direction in isospin space and the subsequent relaxation to the true 
vacuum. We have discussed that the discriminative signature can be provided 
by the two pion correlations in back-to-back momentum configurations ]5]. [T^] . 
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In this paper, we investigate the effect of quantum back reactions onto 
the parametric resonance mechanism. Particle creation, when it occurs due 
to strong coupling with background field oscillations, affects the motion of 
the background fields by acting as dissipation. The back reaction, in turn, 
affects the particle production itself because of the damping of background 
oscillations. We aim at taking account of the interplay between particle pro- 
duction, its back reactions to the background fields, and reaction back to the 
particle production in a self consistent manner. Unfortunately, the effect was 
completely ignored in our previous treatment, but evaluation of its effects is 
indispensable for the correct understanding of a role played by the parametric 
resonance mechanism on DCC formation. This is our first trial toward the 
goal and we explore in this paper the effect of back reactions on single particle 
momentum distributions of pions, with primary concern on its effect on reso- 
nance peaks. We will also calculate correlation length of pions, one of the key 
issues in DCC. 

In our previous works, our formulation relied on the approximation of small 
background sigma oscillations. It may correspond to the final stage of nonequi- 
librium chiral phase transition and we discussed the quantum particle creation 
under the background. We have argued that in this setting we may ignore 
higher order terms in the quantum pion and sigma fluctuations. Due to a 
series of approximations, we are able to write down explicit quantum states 
of pions and sigmas in the form of the squeezed states fL3| . It enables us 
to compute explicitly two pion correlations, and the BCS type back-to-back 
momentum pairing is obtained. It would give a characteristic signature of the 
parametric resonance mechanism. 

In this paper, we calculate the effect of quantum back reactions but still 
within the framework of small amplitude background sigma oscillations. Our 
restriction to the small-amplitude region is not only due to technical reasons, 
but also is physically motivated. To understand the effects of quantum back 
reactions, it is desirable to compare the results with and without back reactions 
at the same initial amplitude of background a oscillation. Of course, the 
limitation clearly implies that the meaning of our results must be carefully 
interpreted, in particular with regard to possible implications to experiments. 

We use the formulation by Boyanovsky et al. M which incorporates the 
quantum back reactions via the Hartree-Fock approximation. In our previous 
computation ignoring back reactions, the pion momentum distributions have 
resonance peaks at the momenta which are characteristic to the parametric 
resonance mechanism. If observed, it would give an unambiguous confirmation 
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of the mechanism. However, it is possible that it can be wiped out when the 
effect of quantum back reaction is taken into account. For related works which 
address DCC by the similar formalism, see e.g., Refs. (14|, [15], |i6f . 

We will observe that the quantum back reactions drastically affect low mo- 
mentum components of the pion field fluctuations but not on the resonance 
peaks. The peak hight is barely changed but the peak position is moved 
to lower momentum. We believe that the fact that the resonance peak sur- 
vives even if the quantum back reactions are taken into account gives us some 
positive hints about experimental detectability of the parametric resonance 
mechanism. 

Our calculation of two pion correlations results in a short correlation length 
of order 1-2 fm which is much shorter compared to the one expected in DCC 
scenario. We will briefly address its possible interpretation, in particular in 
the relationship with the different result obtained by Kaiser [TJj . 

In section 2, we review the Boyanovsky et al.'s formalism. In section 3, we 
calculate the single pion momentum distributions and two pion correlations. 
In the last section, we give concluding remarks. 



2 Quantum Back Reactions in the Hartree- 
Fock Approximation 



We start with the following linear sigma model Lagrangian 

1 

2 

where 



c = -d^ a d^ a -v{ct>), (l) 



V{(f>) = ^(0 a 0a - v 2 ) 2 + ha, a = (a, vf). (2) 
Typical values of the parameters 

A = 20, v = 90 MeV, m w = J— = 140 MeV, (3) 

V 

are employed throughout this paper, the same ones as used in our previous 
analysis. ||. We rely on the formalism based on the Schrodinger picture 
formulation of quantum field theory ||. 

In the following calculations, we consider the pion production in a back- 
ground sigma field which is oscillating along the sigma direction in isospin 
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space. Decomposing sigma model fields into the background and fluctuations, 
a = fo + ip and n = 7? , we substitute them into the Lagrangian (^). We employ 
the Hartree-Fock approximation to take account of the back reaction. Then, 
the quartic terms in the Lagrangian can be replaced by quadratic terms, e.g., 



(p — > 6{cp )(f — 3((p ) . The Hamiltonian is expressed as 
H-- 



+ A(3(^) + (7T 2 )Vo + 



dV{a, 0) 



da 



<T=tfiO I 



(4) 

where II's stand for the conjugate momenta of fields indicated as indices. The 
time dependent effective masses for sigma and pion are written, respectively, 

as 

Ml= AM-^) + 3A(^>, 
Ml= A(^-^) + 3A(7r 2 ). 



(5) 



The expectation value is evaluated by 



( V 2 ) = Tr Up) , 



(6) 



where p is a functional density matrix written in terms of (p. In a self-consistent 
Hartree approximation, the density matrix in sigma sector takes the Gaussian 
form PI 



p a [a, a] = JjA/^exp 



A 



(7)" 

where J\f v , A v , A*, and B v are determined so that the density matrix obeys 
the functional Liouville equation 

i% = [H,p\. (8) 



For instance, we have 



iSf v = lyNyiAy - A 



%k= (3A(^ 2 ) + 3A(7r 2 ))(po + 



tf)l 

dV(a,0) 



da 



AO- 



(9) 
(10) 



°-=<£o / 
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where Q is a finite volume of the system. Actually, the equation ( |T0] ) leads to 
the following equation of motion of y? 

</3 + Mvl - vl)^ + (3A(^ 2 > + 3A(7r 2 ))v9 -h = 0. (11) 

As was done in the analysis by Boyanovsky et al. M, the assumption of 
no cross correlations between sigma and pion fields is adopted in the deriva- 
tion of the above effective Hamiltonian Thanks to this assumption, the 
density matrices of sigma and pion sectors are factorized, allowing their sepa- 
rate treatment. Yet, the interaction between sigma and pion fields influences 
the dynamics through motion of <po, whose behavior is strongly affected by 
quantum back reactions due to particle productions. Among other things, it 
ensures the conservation of total energy despite the factorization of sigma and 
pion sectors. To briefly summarize, the effects of higher order terms of fields 
are included under the Hartree-Fock approximation, and the back reaction 
due to particle production is taken into account through the mean field values 
evaluated with use of the density matrix. 

The expectation value of sigma fluctuation is expressed as 

1 A (12) 

= 22 Tr (<Pk<P-kP*) ■ 
k 

After a short calculation, which is straightforward with the Gaussian density 
matrix, we find 

W- fc ) = ^K| 2 coth(^^) (13) 

under the assumption that the initial state is in thermal equilibrium with 
a temperature T = l//3o- This assumption can be expressed in terms of 
coefficients of the density matrix as follows: 

A v (k,0) = A^(k,0)=w v (k,0)<rth/3ov v (k,0), (14) 
B v (k, 0) = — u v (k, 0)cosech / 9 o; (/ ,(/c, 0), (15) 

AUM) = ( ^^ tanh ^^ ) 1/2 . (16) 



The function ip v defined by 



t—^- = Re.4.^ • tanh /3 u v (k, 0) + ilmA^ (17) 
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obeys the equation of motion of ip, 



i> v (k,t)+u 2 v (k,t)if) v (k,t) = Q, (18) 

with 

u;l(k,t) = k 2 + Ml(t). (19) 
The initial condition for if) v is given by 

Vv(M)= i ) = , ^(fc,0) = i^(fc, 0), (20) 
iC0(p{k, 0) 



through the initial conditions for the density matrix under the assumption of 
thermal equilibrium. The point is that the function ip satisfies the classical 
equation thanks to the transformation (flTf) , despite that it describes the quan- 
tum system. As one can see, the right hand side of the equation ([13]) contains 
the expectation value of if 2 at t = through the effective mass dependence 
in the frequency fliTf). Hence, the initial value of {ip 2 ) must be determined 
self-consistently. 

Let us consider an expectation value of the pion number (n^{k,t)). A 
particle picture in a time dependent background, which corresponds to the 
time-dependent mass @ , has the similar feature as those in theories in curved 
space-time. The effective mass has the similar form as the coefficient of the 
quadratic term of scalar fields coupled with the Ricci scalar ||. In this case, 
the Fock spaces chosen in both of the asymptotic regions (t — > ±oo) belong 
to different Hilbert spaces. A vacuum in one of the Fock spaces is a con- 
densed state of particles defined in the other Fock vacuum. The creation 
and annihilation operators in each Fock space are related via the Bogoliubov 
transformation. Two different definitions of this transformation are employed 
in several literatures so far p|- [|T6[ . 



In general, the Bogoliubov transformation between operators in two Fock 
spaces is written as 

a k (t) = a k (t)a ok + ^(*)<4-fc> 4(*) = Pk(t)ao-k + a* k (t)al k , (21) 

where each of a k and ao k belongs to a different Fock space and the coefficients 
obey \a k {t)\ 2 —\f3 k {t)\ 2 = 1. From the transformation (pT|), the number operator 
a k (t)a,k(t) is of the form 

a\{t)a k {t) = (\a k (t)\ 2 + \(3 k (t)\ 2 ) a\a k + \(3 k (t)\ 2 . (22) 
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A possible choice of coefficients of the Bogoliubov transformation is given by 

(23) 



<*k(t) = I = {tip; + Mk, o)#) ; 

Pk(t) = , (ilPn + Wjf(k, 0)^). 

2vMm) 



and the other possibility is as follows 

a k {t) = 1 (ife + utfct^y J* 

k (t) = 1 {i^ + ^{k ) t)^)e-^ tdt '^ t '\ 



(24) 



The former definition (|23f ) relies on time evolutions of canonical fields, 
e.g., ir(k,t) = U(t)7?(k,0)U~ 1 (t), and Boyanovsky et al. utilizes this type of 
the number operator for the numerical analysis. Here U(t) is a time evolution 
operator. In the latter definition, the creation and annihilation operators diag- 
onalize the Hamiltonian at any time t [17], and a^a defined by these coefficients 
is often referred to as adiabatic number operator. We will choose the second 
option, the transformation coefficients (|4j) as we did in our previous papers. 
By definition, both number operators coincide at t — 0. 

A subtle problem, however, arises in the latter definition. As noticed in 
the previous papers |12|, |J, there exists a possibility that cu(k,t) becomes 
imaginary when the amplitude of the background sigma oscillation exceeds 
a limit, (1 — (po/v) 2 > (m^/m^ 2 . It is highly plausible that the ill-defined 
imaginary frequency is an artifact of the wrong choice of variables, but we do 
not enter into the problem in this paper. 



3 Single Pion Momentum Distributions and 
Two Pion Correlations 

The expectation value of the pion number operator is defined by 

(ntfct)) = Tr(al(t)a k (t)p(0)), (25) 

where p(0) stands for the functional density matrix at t = 0. At t = 0, the 
initial occupation number is given by 

WM)) = ^ ri , (26) 



8 



at finite temperature via the first term of right-hand-side in (^), an expected 
result. It is easy to obtain the energy in pion sector in terms of a function ip^, 
which is defined similarly as ip^ in ([T7|), and it is given by 



E*(k,t) = Q |ife(M)f + lul(k,t) fok(M)| 2 ) I (e?«*M - l) . (27) 

We now proceed to the numerical analysis. We restrict ourselves into small 
initial background oscillation of the sigma field. We parametrize it, throughout 
the remaining sections, by the parameters 

Xo(0) ^ <Po^ (2g) 

V V 

which implies a departure from the bottom of the potential minimum measured 
in units of v, the vacuum expectation value of <p. By solving the coupled 



differential equations (|TT| ) and ([18]), we compute (n^(k,t)). In Fig. 1, we 
present snapshots of the single pion momentum distributions at time t = 0, 
2, 4, 6, 8, and 10 fm with the initial amplitude of a background oscillation 
Xo(0)/f = 0.05. This is the same initial condition as in the previous analysis 
ignoring the back reaction |J and the snapshots obtained there at t = 4, 7, and 
10 fm are presented in Fig. 2 for comparison. (At t = 0, (n^(k,t)) trivially 
vanishes because of the initial condition chosen in ||.) From Fig. 1, one 
observes the three significant characteristics of the time evolution of the pion 
momentum distribution: 

(1) At momenta lower than ~ 100 MeV, the number density grows first and 
then decreases, and seems to undergo damped oscillations. The suppression of 
number density at low momenta seems to indicate that the back reaction due 
to particle creations are efficient at low momenta. It may be a bad news for 
DCC because it may suppress the formation of large coherent domains. We 
will make some comments later. 

(2) At around the resonance energy ~ 200 MeV, a peak starts to develop at 
t ~ 4 fm and continues to grow as time goes by. The series of snapshots in 
Fig. 1 looks like a process that the peak grows by absorbing the ambient pions 
in low momentum region which was originally provided, as seen by comparing 
with Fig. 2, by the background oscillation in an early stage of evolution. 

(3) At momenta larger than ~ 300 MeV, there is virtually no change in the 
initial thermal distribution. 
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The peak height of pion number distribution at resonance without the back 
reaction reaches 0.28 at t = 10 fm in Fig. 2. The corresponding peak height 
with back reaction is 0.38 at t — 10 fm as shown in Fig. 1. The peak height 
with the back reaction is higher but it is mainly due to the fact that the initial 
occupation number (|26|) is nonzero. If we estimate the net increase of height 
due to evolution just by subtracting the initial thermal distribution, we obtain 
0.25. It amounts to about 90% of the peak height without the back reaction. 

We also note, by comparing Fig. 1 with Fig. 2, that the peak position 
moves to lower momentum region when the back reaction is taken into account. 
It can be understood by the following consideration. 

If we ignore back reactions, the equation of motion of pion fields reduces 
to the Mathieu equation 

+ A ~ 2( 1 C0S(2Z)^) 7T fc = 0, (29) 

where A = 4(k 2 + mfyjm % a and q = 4\vxo(fy / m l 0- Here we used a di- 
mensionless time variable defined by z = m a t/2. It is well known that the 
equation admits unstable solutions in a wide range of parameters, in partic- 
ular at A = n 2 for small q which is relevant to our case, where n denotes an 
arbitrary integer [18| . Using the effective pion mass corrected by back reactions 



the peak position is roughly estimated as 



kpeaM ~ \/ ^ - ml ■ ^ - \xl - 3A(tt2), (30) 

for the first resonance band A = 1. 

The back reaction affects k pea k in two opposite ways. The amplitude of the 
oscillation obviously damps by back reactions and it makes the third term in 
the square root smaller in magnitude at t ^ 0, and hence k pea k larger. On 
the other hand, the last term, which was absent in a case of no back reaction, 
tends to let k peak be smaller. A simple estimate shows that the latter effect 
wins. 

To confirm this interpretation, we run the computations with four different 
initial conditions in the region Xo/ v £ [0.04,0.10]. In Fig. 3, we plot the pion 
distributions at t = 20 fm. One notices that the peak position moves toward 
a lower momentum region by increasing the amplitudes of background sigma 
oscillation in consistent with our interpretation. It is also notable that the 
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peak height rapidly increases as Xo/ v § e ts larger. It clearly shows a hint for 
experiments. 

The behavior of damping of the background oscillation is shown in Fig. 
4. This and the resultant increase of the effective mass cause a shift of the 
peak position to a lower momentum, as we just argued. In Fig. 5, we plot the 
time evolution of the energy of pion sector, E n . It is defined as the summation 
over k in fl27|) . It increases while the background sigma oscillation damps with 
characteristic time scale of ~ 10 fm, and becomes stationary at about t ~ 15 
fm when the background sigma oscillations die away. The pion production 
effectively terminates at this point. 

The momentum distribution for sigma fluctuations is shown in Fig. 6. It 
represents a sharp contrast with the behavior of pion distributions, having no 
evolution until t = 30 fm. The first resonance is expected at zero momentum, 
but it is not visible. It may be due to cancellation between the resonance 
enhancement and the damping due to back reactions at low momenta. 

We have observed that the time scale of the energy dissipation is ~ 10 fm. 
The time scale can roughly be understood by the following arguments. For 
simplicity, we consider a system with pion and background fields only. Suppose 
that the background fields oscillate harmonically. Using the virial theorem, the 
energy conservation reads 

m^o(0) 2 + 6,(0) = mlx (z) 2 + e.(0)e 2 ^, (31) 

since the pion field iTk(z) is expressed by using the critical exponent as 7Tfc(0)e^ 
in a resonance band. The pion energy density e n is proportional to the square 
of the pion field 7rfe(z). Xo implies the amplitude of background oscillation, 



and fx is the critical exponent of the Mathieu function [fL8|| . As mentioned 
above, the equation of motion of pion is described by the Mathieu equation in 
this simple system and the energy of pion is dominated by the modes in the 
resonance band in a good approximation. 

In our analysis, a parametric resonance occurs in a narrow resonance band 
and li is given by ii = q/2 |TH], that is 



2\vx 2\v^e -e 7r (z) 
li = — — = § , (32) 

where eo is the initial total energy density of the system considered. Because 
the time dependence of pion energy at any time z is expressed by e n (z) = 
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e 7r (0)e 2/iZ , one obtain the following equation 

de n (z) 



dz 



2^(z) (33) 



in the region where the time derivative of /i can be neglected. Practically, it 
means the narrow resonance region. The solution is given by 



cq I 1 — tanh 2 



2A ^ + ta„lrM-^M 



ml V eo 



(34) 



where z G [0, z en d]- The time z en d is a dimensionless cut off time beyond which 
the equation ([33|) is no longer valid. One can see that almost all energy of 
background field is transferred to pions by the time given by 



m \ -1 / e 7r(0) 

z end = — — —tanh Wl . (35) 

2Xv^M V eo 



The energy density eo is determined, once the initial amplitude of sigma oscilla- 
tion is given. With parameters used in our numerical calculation, Xo/ v = 0.05, 
e — 8.67 x 10 6 MeV 4 . Using the parameters given in (|3]), t en d is estimated 
to be ~ 8 fm. Note that this value is under-estimated because \i is treated as 
a constant in the method we used. Since /i decreases with time, the energy 
transfer needs longer time to be completed. Nevertheless, it is notable that 
the treatment gives a correct order of magnitude estimation of z en d- 

Two point correlation function C(t,r) is similarly obtained by means of 
the density matrix, i.e., 

C(t,r) = (#(t,0)jf(t,r)> = 7>(ff(f,0)jf(f,r)p(0)) 

^E^o(^)Wv| 2 co t h(^)), < 36 > 

where jo(x) = sinx/x denotes the spherical Bessel function of order 0. In Fig. 
7, we plot two point pion correlation functions at t — 5, 10, and 20 fm. At 
t = 10 when the evolution effectively ends, the correlation length £ is given by 
£ ~ 1.3 fm. 

In search for the possibility of correlation length growth, we run the com- 
putations with varying initial conditions. The correlation length increases for 
larger initial background oscillations, but not so much, as we observe in Ta- 
ble. 1. 
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Table 1: Correlation length £ of pions at t = 20 fm for varying initial back- 
ground amplitude (p /v's. 



Xo(0)/v 


0.05 


0.06 


0.07 


0.08 


0.09 


0.10 


£[fm] 


1.27 


1.36 


1.43 


1.50 


1.61 


1.68 



Apparently, this result is in disagreement with the conclusion reached by 
Kaiser who found the large domain formation due to the parametric resonance 



mechanism [[14]]. Our simulation indicates that the back reaction primarily 
affects low momentum components and it is unlikely that the mechanism itself 
leads to large domain formation. 

We note, however, that there are two significant differences between our 
and his calculations: 

(i) The initial amplitude of background sigma oscillation is large, Xo/ v ~ 1 i n 
Kaiser's calculations, whereas it is only up to 0.1 in our present analysis. 

(ii) The quantum back reaction is fully taken into account in each step of evo- 
lution of the system under the Hartree-Fock approximation in our calculation. 
But it is ignored in Kaiser's analysis apart from that the effect is evaluated to 
determine the time at which the simulation has to stop. 

At this stage, we feel it difficult to judge whether or not large coherent 
domains are expected to form by the parametric resonance mechanism. It 
may well be the case that it depends upon the initial conditions. To really 
address the issue, we have to run a simulation with large initial amplitudes 
with full taking account of quantum back reactions. To carry it out, we need 
a formulation which is free from the instability problem as we mentioned at 
the end of Sec. 2. 



4 Concluding Remarks 

In this paper we have investigated features of quantum back reactions due to 
particle production by the parametric resonance mechanism in an environment 
of nonequilibrium chiral phase transition. We have calculated the single pion 
momentum distributions and two pion correlations under the initial small am- 
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plitude sigma oscillation around a potential minimum and with initial thermal 
equilibrium assumption. We relied on the formalism developed by Boyanovsky 
et al. to take account of back reactions under the Hartree-Fock approxima- 
tion based on the Schrodinger picture formulation of quantum field theory. 
Remarkably, we observed that the resonance peak survives under the back 
reaction even in such a strongly coupled linear sigma model. 

If the setting of small initial background oscillation is relevant in physical 
situation in high-energy hadronic collisions, a sharp peak in pion momentum 
distributions can be a clear signal for the parametric resonance mechanism, 
one of the candidate mechanism for DCC 

The potential possibility of large domain formation, the key issue in DCC, 
is not fully explored and hence unanswered in the present analysis. We hope 
that we can return to this problem by having a formalism without instability 
problem in the near future. 
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